Experimental test of strongly non-classical character 
of a noisy squeezed single-photon state 
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We experimentally verify the quantum non-Gaussian character of a conditionally generated noisy 
squeezed single-photon state with positive Wigner function. Employing an optimized witness based 
on probabilities of squeezed vacuum and squeezed single-photon states we prove that the state 
cannot be expressed as a mixture of Gaussian states. In our experiment, the non-Gaussian state 
is generated by conditional subtraction of a single photon from squeezed vacuum state. The state 
is probed with a homodyne detector and the witness is determined by averaging a suitable pattern 
function over the measured homodyne data. Our experimental results are in good agreement with 
a theoretical fit obtained from a simple yet realistic model of the experimental setup. 
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I. INTRODUCTION 

Quantum non-classicality is a fundamental feature of 
quantum physics and a key resource in modern applica- 
tions such as quantum information processing or quan- 
tum metrology. At the early stage of quantum optics, 
coherent states have been identified as quantum states 
establishing a bridge between the classical and quantum 
coherence theory [l|]. Coherent states can be generated 
from vacuum by a classical coherent driving described 
by an interaction Hamiltonian linear in the annihilation 
a and creation operators. Coherent states provide a 
quantum description of classical coherent light waves and 
their mixtures p = J P(a)\a)(a\ d 2 a are therefore con- 
sidered to be classical because they are described by a 
positive semi-definite Glauber-Sudarshan function P(a) 
[l[ , which can be treated as a probability density of wave 
amplitudes. On the other hand, if P(a) cannot be inter- 
preted as a probability density, the state becomes non- 
classical from the point of view of coherence theory. 

Operationally, non-classical states cannot be prepared 
using only coherent states and passive linear optical ele- 
ments. The lowest-order nonlinearity capable of produc- 
ing non-classical state is described by interaction Hamil- 
tonians quadratic in a and Such a quadratic nonlin- 
earity can generate non-classical squeezed states whose 
P(a) cannot be considered as an ordinary probability 
distribution. However, the squeezed states can still be 
described by a positive semi-definite Gaussian Wigner 
function W(a) Q. 

In analogy to the set of mixtures of coherent states we 
can introduce a set Q of all Gaussian states and their 
mixtures, p = J P(X)pa(X) dA, where pc(X) denotes a 
Gaussian state, multi-index A labels all possible Gaus- 
sian states, and P(X) > is a probability density. All 
states in Q possess positive Wigner function due to the 
Gaussian nature. Note that, according to Hudson the- 
orem, Wigner function of any pure non-Gaussian state 



is negative at some points of phase space [2| and these 
states are thus highly non-classical. On the other hand, 
non-Gaussian Wigner function of a mixed state p gener- 
ally does not imply that the state is highly non-classical 
as this non-Gaussianity might arise solely from a non- 
Gaussian distribution, P(X), of states. This is an exam- 
ple of a classical non-Gaussianity. In contrast we define 
quantum non-Gaussian states as all states which do not 
belong to Q [1, 0). The quantum non-Gaussian states 
are highly non-classical because they cannot be prepared 
from thermal or coherent states using only Gaussian op- 
erations such as squeezing, and classical randomization. 
This means that some higher order nonlinearity is neces- 
sarily involved in the state preparation. In view of the no- 
go theorems for entanglement distillation [5H7|, quantum 
error correction Q , and quantum computing [9l-fll| with 
Gaussian states and operations, quantum non-Gaussian 
states therefore represent an enabling resource for contin- 
uous variable quantum information processing [l2l [l3j . 

It is therefore very important to develop reliable tech- 
niques for experimental identification of the quantum 
non-Gaussian states. Determining whether any given 
state is non-classical is generally a very challenging task 
that may require the investigation of an infinite number 
of conditions (l4| . Moreover, many experiments provide 
only partial information about the observed state. For- 
tunately, these difficulties can be overcome by formulat- 
ing specific criteria that provide a sufficient condition for 
non-classicality. Any such criterion can unambiguously 
verify non-classicality of a class of quantum states while 
for other states the result is inconclusive. 

Recently, it has been shown that the quantum non- 
Gaussian character can be witnessed simply by deter- 
mining the probabilities for the occurrence of vacuum 
and single-photon states po — (0|p|0), p\ — 0]- 
This criterion allows us to conclusively certify quantum 
non-Gaussian character of a large set of states including 
those with positive Wigner functions. Simultaneously, 
it can be used to detect the presence of processes with 



2 



higher than quadratic quantum nonlinearity. 

Detection techniques typically register either particle- 
like or wave-like properties of quantum states. In optical 
settings where both types of detectors can be employed, 
a reliable direct measurement of the photon number is 
possible only for a low number of photons. When deter- 
mining the photon number probabilities from coincidence 
measurements with realistic single-photon detectors that 
only distinguish the presence or absence of photons we 
need to ensure that non-classical features are not over- 
estimated 4]. Alternatively, the probabilities po and p\ 
can be estimated indirectly by homodyne tomography 
[l~5l - [l7l ]. A homodyne detector is used to measure ro- 
tated quadratures xe — (ae~ ld +a i e l6 ), where 9 is 
the relative phase between a strong coherent local oscil- 
lator and a signal beam, and the probabilities po and 
Pi are then reconstructed from the measured data by 
post-processing. In addition, coherent displacement or 
even squeezing can be applied to the experimental data 
so that we can determine po and p\ of a displaced and/or 
squeezed version of the original state only by data pro- 
cessing without performing these operations physically 
on the signal mode. This is a very useful technique that 
significantly broadens the applicability of the criterion. 
Highly non-classical character of many states may be 
masked by a Gaussian envelope which can prevent iden- 
tification of quantum non-Gaussian character from the 
knowledge of p$ and p\. The above technique can re- 
move this Gaussian veil and greatly enhance the power 
of the criterion. 

In this paper we experimentally certify the quantum 
non-Gaussianity of a conditionally generated squeezed 
single-photon state [l8l - [25| whose Wigner function is pos- 
itive at the origin of phase space due to noise. The state 
is prepared by conditionally subtracting a single photon 
from squeezed vacuum state and it is probed with a ho- 
modyne detector. Probabilities po and p\ are determined 
by two different estimation methods, namely linear re- 
construction based on pattern functions J26H30J and non- 
linear maximum likelihood estimation |3ll433| . To opti- 
mally witness the quantum non-Gaussian character of the 
measured state we apply a suitable anti-squeezing opera- 
tion to the experimental data. We conclusively prove the 
quantum non-Gaussian character of the prepared state 
with confidence of 2.4 standard deviations. The exper- 
imental results are successfully fitted with a standard 
theoretical model of photon subtraction from squeezed 
vacuum. 

The rest of the paper is organized as follows. The cri- 
terion allowing identification of quantum non-Gaussian 
states is discussed in Sec. II. The experimental setup 
for the generation of photon-subtracted squeezed states 
is described in Sec. Ill and a theoretical model of the 
setup is presented in Sec. IV. The reconstruction of pho- 
ton number distribution from experimental data by pat- 
tern functions and maximum likelihood estimation is dis- 
cussed in Sec. V. Experimental results are presented in 
Sec. VI. Finally, Sec. VII contains brief conclusions. 



II. QUANTUM NON-GAUSSIAN STATES 

A simple and powerful criterion for practical identifica- 
tion of quantum non-Gaussian states has been recently 
proposed in Ref. @. This criterion is based on deter- 
mination of the maximum probability of a single-photon 
state, pi,G) that can be achieved by Gaussian states and 
their mixtures for a fixed probability of vacuum pq. If 
Pi exceeds this bound then the state is quantum non- 
Gaussian. In fact it suffices to maximize p\ over pure 
squeezed coherent states that form extremal points of Q. 
This optimization yields a parametric description of the 
dependence of pi t a on po 
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where r € [0, oo). The boundary curve {T]) is plotted in 
Fig. 1. Points lying above this curve can only be obtained 
for states that do not belong to Q. Interestingly, the 
class of quantum non-Gaussian states is much larger than 
the class of states with negative Wigner function. For 
example, a mixture of vacuum and single-photon states 
with a dominant vacuum contribution, p — p\Q) (0| + (1 — 
p > \, has a positive Wigner function yet it can 
be shown that it cannot be expressed as a mixture of 
Gaussian states for any p > Q . 

Due to the convex structure of the set Q we can con- 
struct a witness of the state's quantum non-Gaussian 
character. The witness is defined as a linear combina- 
tion of po and p\, 



W(a) = ap +pi, 



(2) 



where a < 1 is a parameter specifying the witness [34j . 
The maximum value of W(a) over Q can be found by 




0.4 0.6 

Po 



FIG. 1: The maximum pi achievable by Gaussian states and 
their mixtures is plotted as a function of po (solid blue line). 
The red dashed line represents the non-classicality boundary 
given by Eq. (JS| and the dotted black line indicates the ulti- 
mate physical boundary po + Pi = 1. The dot-dashed black 
line represents a trajectory in (po,pi) plane of a squeezed 
single-photon state (r = 1) subject to varying losses. The 
left-most point corresponds to no losses (j) = 1) while the 
right-most point corresponds to complete losses (rj = 0). 
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inserting formulas (QJ into Eq. ([2]) and maximizing W(a) 
over r. After some algebra one obtains 
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where 



r = -ln- 



a + Va 2 - 10a + 9 



(3) 



(4) 



is the optimal r for a given a. If W^(a) > Wg(<i) then the 
state is quantum non-Gaussian. Each optimal witness 
can be represented by a straight line apo + pi = Wg (a) 
which is tangent to the boundary curve (J]). 

In a similar fashion we can also investigate whether 
the state is non-classical in the sense that it cannot be 
expressed as a mixture of coherent states |a) . In this case 
the maximum p\ achievable for a fixed po is specified 
by Poissonian statistics and no further optimization is 
necessary :35j, 



Pi 



(5) 



where n € [0, oo) is the mean photon number. The 
boundary ([5|) is also plotted in Fig. 1 as a dashed line. 
By analogy, we can use W(a) also as the non-classicality 
witness. The bound achievable by mixtures of coherent 
states reads W c \ = e a ~ x [35|. 

In this work we are interested in the quantum 
non-Gaussianity of approximate squeezed single-photon 
states obtained by photon subtraction from a squeezed 
vacuum [l8- 25]. In the ideal case of perfect single-photon 
subtraction from pure squeezed vacuum we would ob- 
tain pure squeezed single-photon state \ip{ r )) = S(r)\l), 
where the squeezing operation reads 



S(r) 



(6) 



r is the squeezing constant, and x and p denote conjugate 
quadrature operators satisfying canonical commutation 
relations, [x,p] — i. The witness W(a) identifies the pure 
state \tp{ r )) a s quantum non-Gaussian for any amount of 
squeezing, because it is a superposition of odd Fock states 
and po = while p\ > 0. As a more realistic case let us 
next consider a mixed state obtained by sending \ijj(r)) 
through a lossy channel C with transmittance r\. In Fig. 
1 we plot the trajectory of points Po,Pi generated by 
varying the transmittance r\ for a fixed amount of initial 
squeezing r. We can see that the curve enters the region 
of Gaussian mixtures so for high enough losses we can no 
longer identify the state as quantum non-Gaussian or as 
non-classical. There exists a threshold transmittance rjth 
for which our criterion reveals the quantum non-Gaussian 
character of the state. The dependence of 77th on the 
initial squeezing r is plotted in Fig. 2(a). 

Since we deal with a squeezed state, we can intuitively 
expect that a witness based on Fock state probabilities 
will not be optimal. We can significantly improve the per- 
formance of our witness if we first extract the highly non- 
classical core of the state [36[ by anti-squeezing the state, 
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FIG. 2: (a) Threshold transmittances 77th (solid line) and 77th, s 
(dashed line) corresponding to witnesses W(a) and W(a,s), 
respectively, are plotted as function of squeezing constant, (b) 
A (po,Pi) plane trajectory of a lossy squeezed single- photon 
state (r = 0.5, 77 = 0.4) subject to anti-squeezing operation 
with a variable degree of anti-squeezing s (dotted black line). 
The blue solid line and red dashed line have the same meaning 
as in Fig. 1. 



ps = S'(s)pS(s). Note that the squeezing is a Gaussian 
operation that maps mixtures of Gaussian states onto 
mixtures of Gaussian states. Therefore, if ps is quan- 
tum non-Gaussian then p is also quantum non-Gaussian. 
Equivalently, we can introduce a generalized witness 



W(a, s) = ap (s) + pi(s), 



(7) 



where p n (s) — (n\S^ (s)pS(s)\n) are diagonal density ma- 
trix elements in basis of squeezed Fock states S(s)|n). 

The anti-squeezing partly compensates for the initial 
squeezing and makes the state closer to a mixture of 
single-photon and vacuum states for which the witness 
@ is very powerful. Note that since the lossy channel 
and squeezing do not commute, S£(p)S^ 7^ C(SpS^), 
the anti-squeezing of the output mixed state is not fully 
equivalent to reducing the initial squeezing r. In Fig. 
2(b) we plot a typical trajectory of po,Pi pairs when the 
anti-squeezing s is varied. We can see that for a suit- 
ably chosen s we can detect the quantum non-Gaussian 
character of the state although without anti-squeezing 
(corresponding to s — 0) our witness fails. For compar- 
ison we plot in Fig. 2(a) the threshold transmittance 
?7th,s above which the witness W(a, s) detects quantum 
non-Gaussian character of the considered state. We can 
see that 77th, s < %h and this gap is a clear indication 
of the enhanced power and increased applicability of the 
witness W(a,s). Note that in addition to anti-squeezing 
we could also coherently displace the state to remove any 
coherent component. However, this is not necessary for 
our present purpose because the states studied in this 
paper do not contain any coherent component. 



III. EXPERIMENTAL SETUP 

The set-up used for the experimental test of non- 
classical character of photon-subtracted squeezed states 
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|25| is shown in Fig. [3l A cavity-dumped Ti:sapphire 
laser (Tiger-PS, Time-Bandwidth Products) produces 5 
ps long pulses with a maximum energy of 50 nJ, a repeti- 
tion rate of 815 kHz and a central wavelength of 830 nm. 
The laser pulses are up-converted in the process of second 
harmonic generation (SHG) using a 3 mm thick period- 
ically poled crystal of potassium titanyl phosphate (PP- 
KTP1). The remaining radiation at 830 nm is removed 
by a set of dichroic mirrors (DM) . The frequency doubled 
pulses are used as a pump in a similar crystal (PPKTP2), 
phase matched for colincar and fully degenerate paramet- 
ric generation (OPA). It produces a squeezed vacuum 
with squeezing strength tunable from to 3.5 dB. 

The generated squeezed light impinges onto an asym- 
metric beam splitter (ABS) which reflects 7.7% of the 
signal to a single-photon detection setup. It consists of 
a narrowband Fabry-Perot filter (FP, FWHM=0.04nm) 
and a single mode fiber (SMF) that guides the signal 
to an avalanche photo diode (APD) operated in Geiger 
mode (SPCM-AQR-14, Perkin- Elmer) . We estimate the 
total quantum efficiency of the filtering and subsequent 
detection process to be approximately 8±1%. The uncer- 
tainty of the efficiency is quite high particularly due to 
the uncertainty of the detection probability of the APD. 
The signal transmitted by the ABS is mixed with a local 
oscillator (LO) using polarizing beam splitters (PBS1, 
PBS2) and a half- wave plate (HWP), and detected by 
means of a pair of photo diodes (S3883, Hamamatsu). 
The difference of the photo currents is generated and the 
resulting current is amplified by a charge-sensitive am- 
plifier and finally fed into an oscilloscope working in the 
memory-segmentation regime. The acquisition is trig- 
gered by a detection event of APD and a synchronization 
pulse from the laser which suppresses the effect of elec- 
tronic dark counts (effectively below 3 Hz). The homo- 
dyne detector (HD) efficiency of 80±3% is limited mainly 
by the efficiency of the photo diodes itself (94%), the 
mode matching of signal and LO (95%), and the trans- 
mittance of passive optical elements in the signal path 
(95%). 
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FIG. 3: Experimental set-up used to generate photon- 
subtracted squeezed states. 



IV. THEORETICAL MODEL 

To compare the experimental results with theoreti- 
cal predictions we employ a standard model of squeezed 
single-photon state preparation [37T - I39T ] . An equivalent 
scheme of the experimental setup is shown in Fig. 4 and 
our goal is to determine the probabilities pg and p\ as 
functions of model parameters V x , V Pl T, 77, rjn, nth, and 
Q. Our derivation is based on phase-space representation 
and Gaussian-state formalism because the non-Gaussian 
quantum state prepared by photon subtraction can be 
expressed as a weighted difference of two Gaussian states 



The input signal mode A ln is prepared in a generally 
mixed squeezed vacuum state p ln described by a diagonal 
covariance matrix (CM) 7a, in = diag (214, 21^,), where 



V x = ((Ax Am ) 2 ), V p = ((A PA J 2 ) 



(8) 



denote the variances of squeezed and anti-squeezed 
quadratures, respectively, and V X V P > i. A small frac- 
tion of light is then tapped-off by a strongly unbalanced 
beam splitter BS with intensity transmittance T and re- 
flectance R = 1 — T <C 1. We assume that the auxiliary 
input port B of BS is in the vacuum state. The state of 
modes A and B at the output of beam splitter is Gaus- 
sian [!(| and its Wigner function reads, 



1 



it 2 Vdet 7 ab 



(9) 



Here £ = (xa,PA, xb,Pb) is the phase space coordinate 
vector and the two-mode covariance matrix jab is given 
by 



1AB 



Tj Aia + RlB iD VRT('fA ia - 7B in ) 
/ i?T(7A m - 7B in ) RjA ln + Tj Bin 



(10) 




BHD 



FIG. 4: Equivalent theoretical model of the experimental 
setup. V x and V p denote variances of amplitude and phase 
quadratures of the input mixed squeezed vacuum state, T is 
the transmittance of tap-off beam splitter. Inefficient single- 
photon (homodyne) detection is modeled as a sequence of a 
beam splitter with transmittance n (j]h) followed by perfect 
detector. Thermal photons with mean number nth can be 
injected into the signal beam which models electronic noise 
in homodyne detector BHD. Single-photon detector APD can 
be triggered by dark counts or photons not coming from the 
signal mode which occurs with probability 1 — Q. 
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where j Bia = 1 is the covariance matrix of vacuum. It 
is convenient to model inefficient single-photon detection 
as a combination of a lossy channel with transmittance r\ 
followed by a perfect detector with unit efficiency Sim- 
ilarly we can model homodyne detection with efficiency 
r]H and background noise nth by a sequence of a lossy 
channel with added thermal noise followed by a perfect 
homodyne detector. The two-mode covariance matrix 
that accounts for detection inefficiency and noise can be 
expressed as follows, 



Iab = M 1AB M T + G, 



(11) 



where 



M = 





WOO 
^ 
y/rj. 



(12) 



and 



G = 



1 - Vh + 2n th 
















1 — 7/ 





1 



(13) 

The on/off detector placed on mode B can be described 
by projectors onto vacuum and the rest of the Hilbcrt 
space, respectively, n = |0)(0| (no click) and IT = 1 — 
n (click). If a click of the detector occurs the state of 
mode A collapses into the output state 



1 



Pout — 



Tr B (l A ®Tl hB p AB ) 

Tr(l A ®IL hB p AB ) 1-Vo 



(rod) 



(14) 

where p AB is the joint state of modes A and B 
in front of the detectors, p( rcd ) = Tr B (p AB ) is 
the reduced state of the output mode A, Vo = 
Tr (1 A <g) n 0) B Pab) is the probability of no click, and 
= Xr,4 (t A (g) n ,_B Pab) /Vo is the conditional state 
of mode A corresponding to no click of the detector. 

The output state can be most easily calculated using 
the formalism of Wigner functions bacause the POVM 
element Ho is a projector onto vacuum that possesses a 
Gaussian Wigner function. Wigner function of the out- 
put state (fT4")l can be written as a difference of two Gaus- 
sian Wigner functions centered on origin, 



w out (u) 



i 



tt(1 - Vo) 



Vdet 7i 



V dct 7o 



Here £ A — (x A ,p A ) T . In order to express the covariance 
matrices appearing in Eq. (|15j) we split the two-mode 
covariance matrix 7^ B into single-mode blocks, 



The covariance matrices read 7/ = T A and 70 = T A — 
Tc(r_B + l) _1 rj [H and for the probability of no click 
we have, 



Vo 



i) 



(17) 



Despite heavy spatial and spectral filtering, the sin- 
gle photon detector APD can be sometimes triggered by 
photons from other modes than the mode which is subse- 
quently observed by homodyne detector. Such false trig- 
gers can also include dark counts caused by thermal fluc- 
tuations or other effects. We can define an effective mode 
overlap Q as a probability that a click from the APD is 
caused by photons subtracted from the right mode. If the 
APD is triggered by a photon coming from some other 
mode then the output state in mode A is prepared in a 
Gaussian state p' rcd ) with zero mean and covariance ma- 
trix 7/. The overall state in mode A is then a mixture of 
state (fl"4]l with probability Q and a Gaussian state ( o( rod ) 
with probability 1 — Q. The Wigner function of the re- 
sulting state preserves the form (TT5j) . only Vo is replaced 
with 



V'o 



QVo 



l-Po(i-Q)' 



(18) 



Anti-squeezing operation on the output state trans- 
forms the covariance matrix of each constituent Gaus- 
sian component according to jj — > SjjS T , where S = 
diag(e _s , e s ). 

In order to evaluate the witness W(a) we need to calcu- 
late the probabilities po and pi of finding the output state 
(Til)!) in the vacuum state and single-photon Fock state, 
respectively. The probabilities can be derived from the 
overlap formula 



Pk 



2tt 



W out (U)W lk) (U)d 2 U, (19) 



where 



w lk) (U) 



i[2fc£la + (-l) fe ] 



exp ( 



(20) 



are Wigner functions of vacuum (k = 0) and single- 
photon Fock state (k = 1), respectively. Performing inte- 
gration in Eq. (|19p we arrive at the probabilities po and 
Pi in the following form, 



Po 



Pi = 



l-V 
2 



1 



V'o 



v'det (7/ + 1) 
det 7/ — 1 



V'det (70 + 1) 
V'o (det 70-I 



l-V 



[ [det (7/ + 1)] 2 [det (70 + 1)] 



Iab 



r£ r B 



(16) 



In Sec. VI we will use these formulas to find the best 
theoretical fit to the experimental data. 
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V. WITNESS ESTIMATION 

In this section we describe the estimation techniques 
that were used for determination of the probabilities 
Po(s) and Pi(s) from the experimental data. As a main 
tool we utilize the pattern functions that provide unbi- 
ased linear estimators of pj (s) and allow for straightfor- 
ward estimation of error bars, which is particularly im- 
portant in the present case. For the sake of comparison 
we also perform maximum likelihood estimation of p„(s). 



A. Pattern functions 

A homodyne detector measures the probability distri- 
bution w(xg;9) of a rotated quadrature operator xe — 
x cos 6+p sin 9 where 9 is the relative phase between local 
oscillator and signal beam H3] ■ If the measurement 
is performed for values of 9 spanning the whole inter- 
val 9 £ [0, 7r] then the measurement is tomographically 
complete and any element of the density matrix p can be 
determined from the recorded data. We are interested in 
estimation of diagonal density matrix elements in Fock 
basis, i.e. photon number probabilities po and p\. A 
particularly straightforward method is to determine esti- 
mates of p n by averaging appropriate pattern functions 
fn(xe) over the sampled quadrature statistics [26l - l30| . 

1 

Pn = ~ 
IT 



shift 9k — j^ir, k = 1, ...,K. In the actual experi- 
ment, the phase is approximately linearly modulated in 
time and is then fitted and binned into K = 40 discrete 
equidistant values in the [0,7r] interval, so the following 
treatment is applicable to our data. Let Mk denote the 
number of quadrature measurements for phase 9k and 
let Xk, m denote measurement outcomes for this setting, 
m = l,...,Mfc. In experimental data processing, the 
integral (|2ip is replaced with a finite sum, 



1 K 1 

-Y — 



k=l K m=l 



(23) 



The inner sum represents averaging of the pattern func- 
tion f n over sampled quadrature statistics for a fixed 9k , 
while the outer sum represents averaging over the phase 
9. If the number of quadrature samples Mk is not a 
constant then measurement results obtained for different 
values of 9 have different weight in Eq. ([2"B")) . 

Statistical uncertainty of estimated p n can be quanti- 
fied by its variance V(p n ) — (p„) — (p n ) 2 - Since Xk lT n are 
independent random variables, we find that 



1 K 



(24) 



fe=i 



where 



w(x e ;9)f n (xo)dx e d9. (21) V k (J n ) = ^ ^(f&X^)) - (f n (X k>m ))*). (25) 



Note that the pattern functions /„ do not depend on the 
phase 9. For vacuum and single-photon probabilities we 
explicitly have [27i429| 



fo(x) 

h{x) 



2 — 2 v / 7rae x erfi(a;), 

2(2a; 2 - 1) + 4^(1 - x 2 )e~ x2 erti(x), 



™ m— 1 

All quadrature measurement outcomes Xk, m obtained for 
a fixed k exhibit the same statistical distribution, there- 
fore the statistical averages appearing in (j2"5j) can be es- 
timated from the measured data and we obtain 



(22) 



1 M " 1 

V k(fn) = J-p ^2 fn( X k,m) ~ 



™ m— 1 



fn(Xk,m) 



where erfi(x) denotes error function of imaginary argu- 
ment. The functions (|2"21 are plotted in Fig. [5] and we 
can see that they are bounded and asymptotically ap- 
proach in the limit of large \x\. The pattern function 
for estimation of witness W(a) can be obtained as an ap- 
propriate linear combination of the two functions (|22p. 
fw = af + fi . 

Direct reconstruction based on pattern functions is an 
appealing option in situations where we need to deter- 
mine just some particular property of the state such as 
the witness W(a) and do not want to reconstruct the 
whole state. In comparison with more sophisticated sta- 
tistical reconstruction methods such as maximum like- 
lihood estimation, the linear estimation is much faster 
and the statistical uncertainty of any estimated quan- 
tity can be easily determined [4jJ without calculating the 
Fisher matrix or performing complicated Monte Carlo 
simulations. To see this let us assume that the quadra- 
tures were measured for K different values of the phase 



(26) 

It follows from Eqs. (|24| and (|25|) that the variance 
V(p n ) scales as 1/A, where N is the total number of 
measurements. 



2 






1 


' A 
' / 1 
' / I 
i l \ 


\ \ 
\ \ 

— r — \ 1 


~ '~~~-~ r ~ / / \ 




1 






2 







54321012345 

X 

FIG. 5: Pattern functions fo(x) (solid line) and fi(x) (dashed 
line). 
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B. Data anti-squeezing 



Maximum likelihood estimation 



As we have seen in Sec. II, the efficiency of our witness 
can be greatly increased if we construct the witness from 
probabilities of squeezed Fock states. Equivalently, this 
means that we should anti-squeeze the state before we 
estimate the probabilities po an d Pi ■ In the Heisenberg 
picture the anti-squeezing transformation boils down to 
the linear re-scaling of quadratures, 



x e 



P = Poe 



(27) 



This suggests that it may be possible to perform the anti- 
squeezing on the homodyne data without altering the 
experimental setup. In order to show this we express the 
measured quadrature operator xg = x cos 9 + p sin 9 in 
terms of the anti-squeezed quadratures :ro an d Po , 

xg = xoe s cos 9 + poe~ s sin#. (28) 

We can rewrite this expression as follows, 

xe = g(xo cost? + po sin #) = gx&, (29) 

where the new effective phase t? and scaling factor g are 
given by 



tanz? 



-2s 



tan 9, 



(30) 



For comparison, we also perform maximum likelihood 
estimation [3lT - l3"3j of the photon number distribution p n . 
For this purpose we construct a quadrature histogram 
corresponding to averaged quadrature statistics 



w(x) 



1 



c;9)d9. 



(35) 



We divide the real axis into equidistant bins with width 
Aa;. To each measured value Xk t m we assign an index j of 
the corresponding quadrature bin, j = round(X/ Cjrn /Ax), 
and increase the counter for this bin by 1/Mfc, Cj — > 
Cj + 1/Mfc, where initially Cj =0. The resulting val- 
ues Cj are proportional to the probabilities that the 
quadrature outcome falls into the jth window, Pj — 

f^ + i\^ X w(x) dx. We can associate a POVM element IL 

with each bin, Pj — Tr[H jp\. Due to the phase averag- 
ing the POVM elements are diagonal in Fock state basis, 
n j = Er=o n ^»H > 0, and Ej-n.-n = 1. The 
probability of obtaining outcome in jth bin is given by 
Pj = J2 n U-j,nPn an d the likelihood function reads 



L =n e^j 



(36) 



<J 



(31) 



According to the above formulas, the measurement of 
xg can be equivalently interpreted as measurement of 
the quadrature x# of the anti-squeezed state, where 
x$ = xg/g. Photon number distribution p n (s) of the 
anti-squeezed state can be inferred by a modified formula 



Pn{s) = ~ 



The factor 



\w{x e ;6)u{ — ) dx d9. (32) 



,9 Z 



d9 



(33) 



attributes different weights to the pattern function ac- 
cording to the value of 9 because a homogeneous sam- 
pling over 9 is equivalent to an inhomogeneous sampling 
over $ due to the nontrivial dependence of $ on 9. Look- 
ing more carefully at Eq. (|3"2")l we can identify the ef- 
fective pattern functions for reconstruction of diagonal 
density matrix elements in the basis of squeezed Fock 
states, 



f n (xg,0] s) = \f n 

9 



(34) 



Note that these generalized pattern functions explicitly 
depend on 9 through g. 



The maximum likelihood estimates p n which maximize 
the likelihood function L can be numerically determined 
by an expectation- maximization algorithm (32l . |42| . [43| . 

We can generalize the above procedure to estimation 
of p n of the anti-squeezed states. According to the for- 
mula (|29[) . we need to properly re-scale each measure- 
ment outcome, hence the bin index is now given by 
j = round[Afc, m /(gAx)]. Also we must take into ac- 
count the weight factor 4^ = g~ 2 when building the 



phase-averaged histogram, Cj — > Cj + l/{g Mk). In 
this way we obtain a quadrature histogram that corre- 
sponds to the phase-averaged quadrature statistics of the 
anti-squeezed state. Once the histogram is calculated the 
probabilities p n {s) can be determined by the expectation- 
maximization algorithm similarly as before. 



VI. RESULTS AND DISCUSSION 

We have estimated the probabilities po(s) and pi(s) 
from the experimental data using the pattern functions 
described in the previous section. The results are shown 
in Fig.[6l where the probability pairs po(s),Pi(s) are de- 
picted for s e [0,0.4]. The graph suggests that the state 
is both nonclassical and quantum non-Gaussian. How- 
ever, the statistical errors are significant due to a rela- 
tively low number of quadrature samples N = 8000 which 
leads to standard deviations of the order of 1/y/N w 
0.011. 
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FIG. 6: Probability pairs po(s) and pi(s) estimated from ex- 
perimental data with the use of pattern functions (green cir- 
cles). Error bars indicate statistical errors (one standard de- 
viation). Blue solid line represents the boundary above which 
the state is recognized as quantum non-Gaussian. Red dashed 
line indicates the non-classicality boundary and dotted black 
line is the ultimate physical boundary po +pi = 1. 



Since there could be statistical correlations between es- 
timates of po(s) and pi(s), a proper evaluation of the sta- 
tistical significance of the observed non-Gaussian charac- 
ter requires determination of statistical error of the wit- 
ness W(a, s). This is accomplished by calculating the sta- 
tistical error for the sampling function fw as discussed 
in Sec. V.A. We introduce a normalized relative witness 



W R (a, 



W{a,s) - W G {a) 
AW{a,s) : 



(37) 



where AW(a, s) is the standard deviation of W(a, s). For 
each anti-squeezing s we maximize Wr over a. The re- 
sulting optimal witnesses are plotted in Fig. [7[a) and 
the dependence of the optimal value of a on s is shown 
in Fig. [TJb). Without anti-squeezing the confirmation 
of non-Gaussian character is not statistically significant 
but for optimal value of anti-squeezing the quantum non- 
Gaussian character is confirmed with much higher confi- 
dence. The optimum squeezing yielding maximum rela- 
tive witness is s pt = 0.15 and we have 

VF(a opt , aopt) - W G {a opt ) = 0.024 ± 0.010, 

hence the Gaussian threshold is exceeded by 2.4 standard 
deviations. Since po > |, the Wigner function of the 
state is positive at the origin, where we expect it to be 
most negative for this kind of state. The strongly non- 
classical character of the state indicated by the witness 
thus could not be uncovered by looking at the negativity 
of Wigner function. 

Next we compare the experimental results with the 
theoretical model developed in Sec. IV. We optimized 
the model parameters to obtain the best fit to the ex- 
perimental data shown in Fig. [6] In the numerical opti- 
mization we fixed the transmittance of the tap-off beam 
splitter, R = 0.077, and the homodyne detection effi- 
ciency r/H = 0.80, which were both reliably determined 




a 0.74 
0.72 



FIG. 7: (a) The optimal witness W(a op t,s) — Wis(<iopt) and 
(b) the optimal witness parameter a opt are plotted as func- 
tions of the anti-squeezing parameter s. The error bars rep- 
resent one standard deviation. 



by independent measurements. We also fixed the single- 
photon efficiency r\ — 0.08. Although r\ is estimated with 
a large uncertainty, its value mainly influences the success 
probability of the experiment while the shape of the the- 
oretical fit remains practically unchanged when r\ varies 
within the uncertainty interval. The fitting yields the fol- 
lowing parameters: V x = 0.364, V p = 0.705, Q = 0.625, 
and nth = 0. This theoretical fit is compared with exper- 
imental data in Fig. [5] and we achieve very good agree- 
ment between theory and experiment. For comparison 
we also plot in Fig. [8] the results of maximum likelihood 
estimation of po(s) and pi(s). The ML estimates qual- 
itatively agree with the results of linear reconstruction. 
The maximum likelihood method provides slightly higher 
estimates of p± and slightly lower estimates of po than 
linear reconstruction but the difference is less than one 
standard deviation, c.f. Fig. |SJb,c). 



VII. CONCLUSIONS 

In summary, we have experimentally demonstrated the 
quantum non-Gaussian character of a noisy single-photon 
squeezed state generated by conditional photon subtrac- 
tion from pulsed squeezed vacuum. Our approach based 
on a witness constructed from probabilities of squeezed 
vacuum and single-photon states allowed us to reveal the 
quantum non-Gaussian character of the state with rea- 
sonably high statistical confidence. This was achieved 
by an optimization of the variable squeezing parameter 
s. The generalized witness W(a, s) can identify a much 
broader class of quantum non-Gaussian states than the 
original witness W(a) based on the photon number prob- 
abilities po and pi 3,4] because optimization over vari- 
ous Gaussian operations allows one to suitably match the 
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Po s s 

FIG. 8: Theoretical fit (solid black line) to the probabilities po{s) and pi(s) estimated from experimental data using pattern 
functions (green circles). Panel (a) depicts the probability pairs in the po,Pi plane while panels (b) and (c) show the dependence 
of po(s) and pi (s) on s. The grey areas in panels (b,c) indicate one standard deviation of estimates obtained by pattern functions. 
The red triangles represent results of maximum likelihood estimation of po(s) and Pi(s) from the experimental data. The blue 
dashed line in panel (a) indicates the maximum p\ achievable by Gaussian states and their mixtures for a given po while the 
dotted black line marks the ultimate physical boundary po + pi — 1. 



witness to a given state. We emphasize that the witness 
is able to detect the quantum non- Gaussian character of 
states with positive Wigner function. 

We have shown that the pattern functions represent 
an efficient method of estimation of the witness W(a, s) 
from the experimental homodyne data. Moreover, this 
approach provides reliable easy-to-calculate estimates of 
the statistical error bars. The probability pi(s) is an 
overlap of the measured state p with squeezed single- 
photon state. Since this latter state contains only odd 
Fock states in its Fock-state expansion, the probability 
Pi(s) provides a lower bound on the negativity of the 
Wigner function at the origin of phase space. In partic- 
ular, if pi(s) > 0.5 for any s, then the Wigner function 



is negative at the origin. The pattern functions derived 
in the present work thus also provide a useful tool for 
probing the negativity of the Wigner function without 
the need for a full reconstruction of photon number dis- 
tribution. A detailed analysis of this will be the subject 
of a future work. 
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